EuropaJschea Patentamt 

@ Eur °P ean Patont Office 

Office europeen das brevets 




(n) Publication number : 0 61 1 054 A2 



EUROPEAN PATENT APPLICATION 



© Application number : 94300473.9 
(g) Date of filing: 21.01.94 



© int. a.": H03M 13/00 



(So) Priority : 22.01.93 JP 9333/93 
22.01.93 JP 9334/93 

© Date of publication of application : 
17.08.94 Bulletin 94/33 

(S3) Designated Contracting States : 
DE ES FR GB IT NL 

@ Applicant : CANON KABUSHIKI KA1SHA 
30-2, 3-chome, Shimomarako, 
Ohta-ku 
Tokyo (JP) 

© Inventor : Iwamura, Kelichl, c/o Canon 
Kabushlkl Kalsha 
30-2, 3-chome Shimomaruko 
Ohta-ku, Tokyo (JP) 



@ Representative : Beresford, Keith Denis Lewis 
et al 

BERESFORD & Co. 
2-5 Warwick Court 
High Holbom 
London WC1R 5DJ (GB) 



(54) PotynomlnaJ-set deriving apparatus and method. 



I A set of minimal polynomials for generating a 
multidimensional array for decoding algebraic 
geometric codes is derived at a high speed. In 
order to obtain a set of minimal polynomials F 
for generating a given multidimensional array, 
when sequentially updating a set of polynomials 
F, df n M are not directly calculated, and a newly 
introduced set of polynomials B and the set of 
polynomials F are updated using the highest- 
degree coefficient d, of polynomials which be- 
long to the set of polynomials B. An array 
memory for storing a given multidimensional 
array u. and first and second polynomial 
memories for storing the set of polynomials F to 
be obtained and a set of auxiliary polynomials 
Q, respectively, are provided. In the calculation 
of polynomials f» and df n+1 <*> , accessing oper- 
ations tor respective memories and accessed 
addresses are controlled in parallel depending 
on the degrees of polynomials fl*. 
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This invention relates to error correction, in which in a digital communication system and a digital storage 
system, errors produced in a communication channel or a storage medium are corrected at the reception side 
using error correction codes. More particularly, the invention relates to a polynomial- set deriving apparatus 
and method, in which in a decoding operation using algebraic geometric codes as error correction codes, a 
5 set of minimal polynomials for generating a given multidimensional array corresponding to a syndrome in re- 
ceived words is obtained. 

RS (Reed-Solomon codes) codes and BCH (Bose-Chaudhuri-Hocquenghem) codes are well known as er- 
ror correction codes for correcting errors produced in a communication channel or a storage medium at the 
reception side in a digital communication system and a digital storage system. These codes are actually used 
10 In apparatuses which deal with compact discs, satellite communication and the like. 

Recently, however, codes named algebraic geometric codes utilizing the algebraic curve theory have been 
widely studied (see references (1) - (3)). Algebraic geometric codes belong to a code system which has a very 
wide applicational range, and which includes the above-described RS codes and BCH codes. It has gradually 
become clear that this system includes new codes having better properties than those of conventional codes 
is (see reference (4)). 

Decoding methods have always been a problem while such excellent codes have been discovered, and 
development of an efficient decoding algorithm has been requested. However, an efficient decoding algorithm 
has not been otained at all until Justesen et al. have proposed a generalized Peterson algorithm (see reference 
(5)) in 1989. Justesen et al. have shown that by applying the Sakata algorithm (see reference (6)) proposed 
20 by Sakata to their generalized Peterson algorithm to provide an efficient algorithm for obtaining a two-dimen- 
sional linear feedback shift register, having a minimum number of storage elements, for generating a two-di- 
mensional array having a given finite size, an error locator function can be derived at a high speed with a small 
amount of calculation. 

However, since the generalized Peterson algorithm has a limitation in its correction capability, Skorobo- 
28 gatov et al. have proposed a Modified decoding algorithm (see reference (7)) which guarantees a higher cor- 
rection capability. 

Furthermore, since the Sakata algorithm includes a large amount of useless calculation not directly relating 
to the Modified decoding algorithm, Kamiya et al. have proposed an algorithm obtained by modifying the Sa- 
kata algorithm so as to conform to the Modified decoding algorithm in 1992 (see references (8) and (9)). 

30 On the other hand, the Berlekamp-Massay algorithm (BM method) and the Euclidean algorithm (Eu meth- 

od) are well known as conventional decoding methods for RS codes and BCH codes. The above-described 
generalized Peterson algorithm is an extension of the Peterson method used for decoding RS codes. The Sa- 
kata algorithm is also called a two-dimensional BM method, and is well known as an extension of the BM meth- 
od used for decoding RS codes (hereinafter termed a one-dimensional BM method). 

35 However, an algorithm corresponding to an extension of the Eu method used for decoding RS codes (here- 

inafter termed a one-dimensional EU method) has not previously been considered. 

Although Sakata has also proposed a multidimensional BM method (see reference (10)) as an extension 
of the two-dimensional BM method, a multidimensional Eu method corresponding to that method has not been 
devised. 

40 The multidimensional BM method is expressed by an algorithm having a structure as shown in FIG. 2 (see 

references (6) and (10) for a method of determining a defining point, rs, rt and the like). 

Point n shown in FIG. 2 is updated by ordering termed a total order (see references (6) and (10)). In FIG. 

2, the calculation of df n fl> in step S22 and the updating of f„« using h in step S25 are sequentially performed. 

This is because the calculation of df n o in step S22 at point n is performed using h obtained in step S24 at point 
45 n - 1, and the updating of f n (D using h in step S24 at point n is calulated using df„«> in step S22 at point n. In 

this case, as shown in FIG. 12, a useless time is produced in the processing time, and therefore such a method 

is not efficient 

On the other hand, a method has been proposed, in which the processes of steps S22 and S24 are exe- 
cuted in parallel utilizing that processing for one-dimensional variables in the one-dimensional BM method, 

so which Is a decoding method for RS codes and BCH codes, is sequentially performed (see reference (1 1 )). How- 
ever, since this method uses shift registers having a f ixed number of stages, (t+1 ) (t is the maximum value of 
the defining point at the final point (n ■ p: see FIG. 2) of polynomial f n W) processing clock pulses are always 
required for a single updating operation (cycle) of f n w. 

Accordingly, the apparatus of reference (9) requires useless processing dock pulses in the mid course of 

55 calculation (n < p) of f n f> which is not the largest defining point Furthermore, the method of reference (9) does 
not consider the multidimensional BM method, and therefore is unsuitable for the multidimensional BM method 
in the following points: 

1) Only one polynomial belongs to each of sets of polynomials Fand G in the one-dimensional BM method, 
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but a plurality of polynomials belong to each of sets of polynomials F and G in the multi-dimensional BM 
method. 

2) In the one-dimensional BM method, since the degrees of polynomials which belong to sets of polyno- 
mials F and G and the arrangement of a given array are one dimensional (one variable), respective coef- 

5 ficients can be sequentially stored in a one-dimensional memory (a shift register or the like) from the high- 

est degree. However, since the degrees of polynomials and the arrangement of the array in the multidi- 
mensional BM method are not one-dimensional, efficient processing cannot be performed with a one-di- 
mensional memory. 

3) In the one-dimensional BM method, since only one polynomial belongs to each of sets of polynomials 
10 F and G, parallel processing of polynomials is meaningless, in the multidimensional BM method, however. 

since a plurality of polynomials belong to each of sets of polynomials F and G, parallel processing of cal- 
culation of polynomials has a meaning. 

4) In the one-dimensional BM method, since only one variable is present in a polynomial which belongs 
to each of sets of polynomials F and G, the processes of steps S22 and S23 cannot be performed in parallel 

15 if calculation for one variable is performed in parallel. In the multidimensional BM method, however, since 

a plurality of variables are present in a polynomial which belongs to each of sets of polynomials F and G, 
the processes of steps S22 and S24 can be performed in parallel if calculation for one variable is performed 
in parallel. 
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It is an object of the present invention to shorten the overall processing time by executing a plurality of 
processing in parallel when deriving a set of minimal polynomials forgenerating a given multidimensional array. 

45 Different aspects of the present invention provide a two-dimensional Eu method corresponding to an ex- 

tension of the one-dimensional Eu method, and provide a multidimensional Eu method corresponding to the 
multidimensional BM method by extending this two-dimensional Eu method. 

In certain aspects, the invention provides a method for efficiently performing parallel processing of a plur- 
ality of processes in the multidimensional BM method, and an apparatus for realizing the method. 

so One embodiment of the present invention relates to a polynomial-set deriving apparatus for obtaining a 

set of minimal polynomials forgenerating a given multidimensional array, comprising first storage means for 
storing a first set of polynomials to be obtained, second storage means for storing a first set of aixi liary poly- 
nomial for the first set of polynomials, third storage means for storing a second set of polynomials different 
from the first set of polynomials, fourth storage means for storing a second set of auxiliary polynomials for the 

55 second set of polynomials, first discrimination means for discriminating if the coefficient of a predetermined 
degree of each polynomial of the second set of polynomials stored in the third storage means equals zero, de- 
termination means for newly determinating a defining point if a polynomial, in which the coefficient of a pre- 
determined degree does not equal zero, is present as a result of the discrimination by the first discrimination 
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means, first derivation means for deriving polynomials which belong to the second set of polynomials based 
on the value of the defining point determined by the determination means, the second set of polynomials stored 
in the third storage means, and the second set of auxiliary polynomials stored in the fourth storage means, 
first updating means for erasing all polynomials, in which the coefficient of a predetermined degree does not 

s equal zero, from the third storage means, and for storing the polynomials derived by the first derivation means 
in the third storage means, second derivation means for deriving polynomials which belong to the first set of 
polynomials based on the value of the defining point determined by the determination means, the first set of 
polynomials stored in the first storage means, the first set of auxiliary polynomials stored in the second storage 
means, and the coefficient of a polynomial of the second set of polynomials relating to the determined defining 

10 point, second updating means for erasing all polynomials corresponding to the polynomials erased by the first 
updating means from the first storage means, and for storing the polynomials derived by the second derivation 
means in the first storage means, second discrimination means for discriminating the presence of a change 
of the defining point, third updating means for updating the first set of auxiliary polynomials stored in said sec- 
ond storage means based on the first set of auxiliary polynomials stored in the second storage means, and 

is the polynomials erased by the second updating means, when the second discrimation means discrimated the 
presence of a change in the defining point and fourth updating means for updating the first set of auxiliary 
polynomials stored in the second storage means based on the second set of auxiliary polynomials stored in 
the fourth storage means, and the polynomials erased by the first updating means, when the second discrim- 
ation means has discrimated the presence of a change in the defining point 

20 Another embodiment relates to a method for obtaining a set of minima! polynomials for generating a given 

multidimensional array, comprising the steps of setting initial values for a first memory for storing a first set 
of polynomials to be obtained, a second memory for storing a first set of auxiliary polynomials for the first set 
of polynomials, a third memory for storing a second set of polynomials different from the first set of polynomials, 
and a fourth memory for storing a second set of auxiliary polynomials for the second set of polynomials, dis- 

25 criminating if the coefficient of a predetermined degree of each polynomial of the second set of polynomials 
stored in the third memory equals zero, newly determining a defining point if a polynomial, in which the coef- 
ficient of a predetermined degree does not equal zero, is present as a result of the discrimination, deriving poly- 
nomials which belong to the second set of polynomials based on the value of the determined defining point 
the second set of polynomials stored in the third memory, and the second set of auxiliary polynomials stored 

30 In the fourth memory, erasing all polynomials in which the coefficient of a predetermined degree does not equal 
zero, and updating the third memory by storing the polynomials derived by first derivation means, deriving poly- 
nomials which belong to the first set of polynomials based on the value of the determined defining point the 
first set of polynomials stored in the first memory, the first set of auxiliary polynomials stored in the second 
memory, and the coefficient of a polynomial of the second set of polynomials relating to the determined defining 

35 point erasing all polynomials corresponding to the polynomials erased by the updating operation of the third 
memory, and updating the first memory by storing the derived polynomials, discriminating the presence of a 
change of the defining point, updating the first set of auxiliary polynomials stored in the second memory based 
on the first set of auxiliary polynomials stored in the second memory, and the polynomials erased by the up- 
dating operation of the first memory, when it has been discriminated that the defining point has changed by 

40 the discrimination operation, and updating the first set of auxiliary polynomials stored in the second memory 
based on the second set of auxiliary polynomials stored in the fourth memory and the polynomials erased by 
the updating operation of the third memory, when it has been discriminated that the defining point has changed. 

Still another embodiment relates to a polynomial-set deriving apparatus for obtaining a set of minimal poly- 
nomials for generating a given multidimensional array, comprising array storage means for storing a given mul- 

45 tidimensional array u, first polynomial storage means for storing a set of polynomials F to be obtained, second 
polynomial storage means for storing a set of auxiliary polynomials G for the set of polynomials F, first cal emu- 
lation means for obtaining polynomials f» which belong to the set of polynomials F based on polynomials fro 
stored in the first polynomial storage means, polynomials g® stored in the second polynomial storage means, 
and deviations of polynomials df n n> , second calculation means for obtaining deviations of polynomials df^oo 

so based on the coefficients of the polynomials fl*> calculated by the first calculation means and the multidimen- 
sional array u stored in the array storage means, and control means for controlling accessing operations for 
the first polynomial storage means and the second polynomial storage means and accessed addresses in par- 
allel depending on the degrees of the polynomials fw. The calculation by the first calculation means and the 
calculation by the second calculation means are executed in parallel. 

55 Yet another embodiment relates to a method for obtaining a set of minimal polynomials for generating a 

given multidimensional array, comprising the steps of storing a given multidimensional array u in an array mem- 
ory, obtaining polynomials fw which belong to a set of polynomial F to be obtained, based on polynomials fw 
stored In a first polynomial memory for storing the set of polynomials F, polynomials go stored in a second poly- 
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nomial memory for storing a set of polynomials G different from the set of polynomials F. and deviations of 
polynomials df„». obtaining deviations of polynomials df„^ » based on the polynomials f»> obtained by the first 
calculation operation and the multidimensional array u stored in the array memory, and controlling accessing 
operations for the first polynomial memory, the second polynomial memory, and accessed addresses in par- 
s allel depending on the degrees of the polynomials iw and executing the first and second calculation operations 
in parallel. 

Other objectives and advantages besides those discussed above shall be apparent to those skilled in the 
art from the description of specific embodiments of the invention which follows. In the description, reference 
is made to the accompanying drawings, which form a part thereof, and which illustrate examples of the inverv 
10 tion. Such examples, however, are not exhaustive of various embodiments of the invention, and therefore ref- 
erence is made to the claims which follow the description for determining the scope of the invention. 

BRIEF DESCRIPTION OF THE DRAWINGS 

is FIG. 1 is a multidimensional-array generation circuit according to the present invention; 

FIG. 2 is a multidimensional-array generation algorithm according to a conventional BM method; 
FIG. 3 is a multidimensional-array generation algorithm according to the Eu method of the present inven- 
tion; 

FIGS. 4(a) and 4(b) are diagrams illustrating the operation of the conventional BM method; 
20 FIGS. 5(a) and 5(b) are diagrams illustrating the operation of the Eu method of the present invention; 

FIG. 6 is an algorithm of a multidimensional BM method according to the present invention; 
FIGS. 7(a) through 7(c) are diagrams illustrating the operation when the algorithm shown in FIG. 6 is exe- 
cuted; 

FIG. 8 illustrates a first embodiment of the algorithm shown In FIG. 6; 
25 FIG. 9 illustrates a second embodiment of the algorithm shown in FIG. 6; 

FIG. 10 is a diagram illustrating address division in a memory; 

FIG. 11 illustrates a third embodiment of the algorithm shown in FIG. 6; 

FIG. 12 is a diagram illustrating the operation when the algorithm shown in FIG, 2 is executed; 

FIG. 13 is a diagram illustrating the configuration of an apparatus for decoding algebraic geometric codes; 
30 and 

FIG. 14 is a flowchart illustrating a decoding procedure of the decoding processing unit shown in FIG. 1 3. 
DESCRIPTION OF THE PREFERRED EMBODIMENTS 

35 FIG. 1 3 is a diagram i llust rating the configuration of an apparatus for decoding algebraic geometric codes. 

In FIG. 13, input unit 131 includes a reception device for receiving data from satellite broadcast or a commu- 
nication network, a reading circuit for reading data from a storage medium, such as a CD (compact disc) or 
the like, and the like, and inputs received words corresponding to image data or voice data. Decoding proc- 
essing unit 1 32 decodes the received words input to input unit 131. Output unit 1 33 outputs decoded data, and 

40 includes a display for displaying image data, a speaker for outputting decoded voice data, and the like. 

FIG. 14 is a flowchart Illustrating a decoding procedure of decoding processing unit 132. In the present 
embodiment, as in the case of decoding RS codes, decoding of algebraic geometric codes is performed by 
the following procedure when erasure correction is not included. First in step S141, a received-word series 
having a predetermined size, serving as a unit of decoding, is input In step S142, a syndrome polynomial u(z) 

45 is generated from the received-word series. In this step, a multidimensional array u corresponding to the coef- 
ficients of the polynomial is generated. In step S143, a set of error locator polynominals F and a set of error 
value polynomials B are derived from the multidimensional array u. In step S144. an error is corrected at the 
position of the error in the received words obtained from the set of error locator polynomials F based on the 
value or the error obtained from the set of error value polynomials B. 

so a method and apparatus for deriving the set of error locator polynomials F and the set of error value poly- 
nomials B from the multidimensional array u wiD now be described. 

For that purpose, a two-dimensional Eu method corresponding to an extension of the one-dimensional Eu 
method will be proposed. 

It will be shown that the two-dimensional Eu method outputs results equivalent to those of the Sakata al- 
55 gorithm, which is the two-dimensional BM method, and the Kamiya algorithm, which is a modified algorithm 
of the Sakata algorithm. Furthermore, a multidimensional Eu method corresponding to the multidimensional 
BM method will be proposed by extending the two-dimensional Eu method. 

It will also be shown that the multidimensional Eu method of the present invention is an algorithm which 
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has effects different from those of the conventional multidimensional BM method in providing an apparatus 
and high speed, and differences between the two methods wiil be described. 

The known conventional algorithm of the multidimensional BM method has the structure shown in RG. 2. 
On the other hand, the algorithm of the multidimensional Eu method of the present invention has the structure 

5 shown in FIG. 3. The multidimensional BM method differs from the multidimensional Eu method in that in the 
multidimensional BM method, df n <» are directly calculated in step S22, and the set of polynomials F is updated 
using the calculated values in step S25, while In the multidimensional Eu method, df n to are not directly calcu- 
lated, and the set of polynomials B is updated in step S34 and the set of polynomials F is updated in step S35 
using the highest-degree coefficient di of polynomials which belong to the newly introduced set of polynomials 

10 B. 

The above-described references (8) and (13) have proved that the above-described conventional multidi- 
mensional BM method devives a set of minimal polynomials F for generating a given miltidimenslonal array 
u. Accordingly, by proving that the multidimensional Eu method of the present invention devives the same set 
of polynomials F as that of the multidimensional BM method by the following theorems 1 through 3, it will be 
is shown that the set of polynomials F obtained by the multidimensional Eu method shown in FIG. 3 is also a set 
of minimal polynomials for generating a given multidimensional array u. The multidimensional Eu method of 
the present invention will now be more specifically described illustrating the known conventional two-dimen- 
sional BM method and multidimensional BM method. 

20 First Embodiment 

First, the Sakata algorithm will be shown after the following preparation (explanation of terms). Since the 
proof that this algorithm derives a set of polynomials F for synthesizing a two-dimensional linear shift register 
having a minumum number of storage elements for generating a given finite two-dimensional array u has been 
25 made in the above-described reference (8), an explantaion of the proof will be omitted. 

Preparation 1 (for details, see reference (8)) 

Z A set of pairs n = (n1, n2) of all non-negative integers n1 and n2. 
30 n : This notation is called a point, which is identified with a point having coordinates (n1, n2) on the X- 

Y plane. 

<T : This notation is called a total order, and determines the the magnitude relationship of point n on 
set £. The point next to point n - (n1 , n2) with respect to <T is defined in the following manner 

35 n + 1 : - (nl - 1, n2 + l) (when nl > 0 ) 

: (n2 + 1, 0) (when nl = 0 ) . 

40 < p : This notation is called a semiorder, which is defined as follows: 

m^pn when and only when ml ^ n1, m2 ^ n2 ( 
where m = (ml, m2) and n = (n1, n2). 
m <p n indicates that m s p n and m * n. 

EtP: = {mcE 1 1 ^ p m, m < T p}. 

45 u: u is a finite partial two-dimensional array having a size q, and is defined as a mapping onto field K 

fromXo*. 

F: Brvariate polynomials on field K are expressed by f = Xfrn z™. 

merf, 

and 

so z m = x ro1 y* rf = {m g £ | f m * 0 } 

8 = LP(f) = max {ml men). 
A set of polynomials is expressed by: 

F = {f*°>, fli * 

df n O : If the mapping from point n on set £ onto field K is represented by un, df n ® = ItfmTOum+n-s for poly- 
65 nomials fTO satisfying LP(fW) = sfl>. 

merf. 

V(u): If the set of un for p St q is expressed by: 

u p = {un n e S^}, and If df„W = 0 (0 s£ n < p) 

6 
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for up, f [up] = 0. At that time, it is assumed that 

V(up) = {f(polynomial)|f[uP] = 0}. 
Defining point: When the set of polynomials F is a set of minimal polynomials which can generate a given 
two-dimensional array. LP(fft ) = s<o is termed a defining point A set of minimal polynomials satisfies the fol- 
lowing conditions: 

0 FgV(up) 

ii) i and j satisfying conditions that 1 £ i, j £ 1 , i * j, and LP(fW) ^> LP(fli>) are not present 
Hi) A polynomial g satisfying conditions that g V(u) and LP(g)eAF is not present where 

AF = r/u£ |3 (k), andZs(k) = { n *Z|a(k) < p B >. 



A: This symbol indicates AF when the above-described conditions are satisfied, and is obtained in the fol- 
lowing manner 

15 A = U As - » 

q, e gZqP, 

and A«r ■ = {m e 2 1 m £ p q-s} when f eV(u<?) satisfying conditions that d^, * 0, and LP(f) = s is present 
A«r» = * in other conditions. 

hi: hi of type <i, j) is defined as follows: 
20 hi = Z<- •C0-f0> - (d/d^- " + ™ • ^-gO), 

where 

r=(rt,r2), 

ii = max{s10>. n1 - s1© + 1} f 

r2 = max{s2( | >, n2 - s20*D + 1}, tfl> = LP(ga)) f fffl G V (u") t g<D G V(u»") t d, = df n fl>, and 
25 dj = dg ro a>. 

O: This notation is called a set of auxiliary polynomials of F, and is a set of polynomials 

G = few, - , gd • *>}. 
Lemma 1: j of the type determined by 0 is j satisfying the condition of LP(flD) = s®, 

andp - sW^pSKD - l) ( a20*D - 1). 
30 Lemma 2: k of the type determined by © is k satisfying the condition of LP(fW) < p t 

Sakata Algorithm 

1) n = (0,0), F = {1},andG=*. 

2) Calculate df n 0) for all polynomials of F. 

3) If f(i) satisfying df n tt * o is present a new A and a new defining point t are determined. 

4) The following procedure is executed for all defining points t 
(D t = (s1fl>, s2M) f Provide polynomial hi of lemma 1 . 
® t = (n1 - a1© + 1, n2 s20*i> + 1) ~» Provide hi of type <k, i> of lemma 2. 
© t = (n1 - a10> ♦ 1, s2<J>), 1 £ I s 1 - 1, -> Provide hi of type <j, i>, 
® t = (ai©, n2 - s2© + 1), 2 ^ j si 1 Provide hi of type <i, j-1>. 
® t = (n1 + 1, s2U>) hi = x" 1 -» 1 <B* 1 .fl> 

© t = (310> ( n2 + 1) hi = yr*-«2(9*1.fCI). 

All polynimials satisfying the condition df n o * 0 from F, and all newly obtained hi are inserted in F. 

5) When A has changed, polynomials of set of auxiliary polynomials G for the new F are selected from 
polynomials of the old G and the polynomials removed from F. v 

6) The process is termined ifn = n + 1;n = p. The process proceeds to 2) in other cases. 
A new algorithm according to the present invention which derives the same output as the Sakata algorithm 

will now be described. 



Algorithm 1 



1) n = (0. 0). F - {1}, G = *, A = and B « {u}. 
5 5 2) If the highest-degree coefficient of all polynomials of 6 does not equal zero, a new A and a new defining 

point t are determined. 

3) The following procedure is performed for ail the defining points t 

® t = (siw, s2W) -» Provide polynomial h2 having the type of lemma 1. 
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© t = (n1 - siw + 1 , n2 - s2^ 1 ) ♦ 1) -» Provide h2 of type <k. l> of lemma 2. 
® t - (n1 - siro + 1. s2«). 1 S I m 1 - 1 Provide h2 of type <j. i>. 
® t = (sin). n2 - s2© ♦ 1). 2 s j 2S 1 Provide h2 of type <i, j-1>. 
® t = (n1 + 1. s2<&) h2 = r <M-»i<D*i>.b© 
5 © t = (3lo > n2-H)->h2=y^-«w- 1 .b<o. 

All polynomials in which the highest-order coefficient does not equal zero are removed from B, and 
. all newly obtained h2 are inserted In B. 

4) Polynomial hi of the type determined in 3) is provided. Polynomials corresponding to the polynomials 
removed from B are removed from F. and newly obtained hi is inserted in F. 
10 5) When A has changed, polynomials of set of auxiliary polynomials G for the new F are selected from 

polynomials of the old G and the polynomials removed from F. Polynomials of set of auxiliary polynomials 
A for the new B are selected from polynomials of the old G and the polynomials removed from B. 
6) The process is terminated ifn=n + 1;n = p. The process proceeds to 2) in other conditions, where A 
and B are sets of polynomials such that A = {a<°>, ••• , a< 1 * 3} and B =» {bw, »• b< 1 - 
is h2: h2 of type <i, J> is defined as follows: 

h2 = z r - »d>-b<i) - (dj/dj)*' ■ n + m - W-afl). 
A large difference between the algorithm 1 and the Sakata algorithm is that in the Sakata algorithm. df„» 
are directly calculated and it is checked If all of them equal zero, while In the algorithm 1 , df n t0 are not directly 
calculated but it is checked if the highest-degree coefficient of polynomials b® which belong to B equals zero. 
20 For that purpose, the algorithm 1 introduces sets A and B of polynomials which are not present in the Sakata 
algorithm. 

The fact that set F obtained by the algorithm 1 is the same set as set F obtained by the Sakata algorithm 
can be proved in the following manner 

25 Theorem 1 

If a polynomial obtained by multiplying polynomial u having the coefficient un (neSOp) defined by expres- 
sion (1) by polynomial f in which LPO » s is represented by b, the f- "degree coefficient bv-rr+s of polynomial 
b is expressed by the following expression (2): 
30 u = Zunzv - n (1) 

ne2o p 

bv - n + s = £fhvum + n - s (2) 

merf, 

where v is an arbitrary integer. 



Proof: 



b = f«u = Z ^ fm-ui-zv-i+m 



40 from expression (1). 

Hence, the coefficient bv-n+s of z^ ■*» of polynomial b becomes as shown in expression (2), because i = m + 
n - s from v-i + m = v- n + a. 

(Proof completed) 

45 

From this theorem, it can be understood that the coefficient of 2T of b{i) - fTO-u corresponding to f (i) e F 
equals df n <9 of the Sakata algorithm. 

If polynomial f(i) which is V(u n ) is represented by f n (9, the following theorem holds: 

so Theorem 2 

If polynomials updated by expressions (3) and (4) are represented by l q « and c q W (q > T n >y m), and poly- 
nomial w q W is defined by expression (5), w q M can be updated by expression (6): 

l q « = z».| n W - (d/dj) ^© (3) 
55 CqW = z~.c n n> - (d^dj) zrt.C m (I> (4) 

w q oo = l q (*>*e - Cq^-d (d and e are arbitrary polynomials) (5) 
w q « = z«<w n O - (d/dj ).z<t.w m CD (6), 
where rs and rt are arbitrary integers, and the initial values of l n W, i m o>, cfl and C m <D are arbitrary polynomials. 
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Proof: 

Since expression (6) holds for arbitrary q and k, the following expressions hold: 

w n <0 = l n «. e - CnO-d 

s w m © = l m <D. e - CJD.d. 

Accordingly, the following expression holds: 

Wq ( k ) = 1 Q < k } . e _ Cq < k > • d 

10 = (z»« -loCi) - (di/dj ). z *t 1»( J) )- e 

- (z rfl -cn(j) - (di/dj )'z** -c-ti) )-d 
= z". Wn (t) - (di/dj J) . 

(Proof completed) 

In the Sakata algorithm, IW = f q 0<) g F which becomes V<u<i) is updated by h 1 , and g(D equals fd> = f m 0> which 
20 becomes V(u m ). Hence, hi can be expressed in the form of expression (3). 
The following theorem also holds: 

Theorem 3 

25 If l q » = y*> ,e = u,d = ■£* x, c.^ = 1 , and c^O = 0 in expression (5), w q « = b q « at deg w q « ^ V. At that 

time. b q <*> is the (v-q+s<*>)-degree polynomial. 

Proof: 

30 Since l q w = f q w and e ■ u. expression (5) becomes w q w = b q w - CqW-d. If deg cl^ s 0 and deg Co^ s 0, 

c q w is a positive-degree polynomial from expression (4), and d = z v x, and therefore deg CqW.d ) v. Hence, if 
deg w q (k) ^ v. w q <*> ■ b q <*>. 

From the definition, the highest degree of b q w = u-f q <*> is v + sfr>. Since f q ro at that time is f q We V(uq), dtp* 
=» 0 (I ■ 0, ~. q-1). From theorem 1, df^o is the coefficient of z*-***). Hence, the (v- I + s<*)-degree (i = 0, . q- 
35 1 ) coefficient of b q w becomes 0. and the highest degree of b q w becomes v - q + s«. 

(Proof completed) 

Accordingly, b q <*> = w q « can be updated by expression (6). That Is, b q w can be updated by h2. In this case, 
40 from theorem 1 and theorem 3, di and dj are the coefficients of z*-"*^ and z*-™*^ of polynomials b 0 <» and 
b m <tt, respectively, that is. the highest-degree coefficients. 

As described above, it was proved that F obtained by the algorithm 1 is the same set of polynomials as F 
obtained by the Sakata algorithm. Since it has been proved in reference (8) that the set of polynomials F ob- 
tained by the Sakata algorithm is a set of mininal polynomials for generating a given two-dimensional array u, 
45 it can be said that the set of polynomials F obtained from the algorithm 1 is also a set of minima) polynomials 
for generating a given two-dimensional array u. 

Accordingly, this algorithm can be realized, for example, by the apparatus shown in FIG. 1. First, it is de- 
termined if the highest-degree coefficient of polynomials B stored In memory 13 for storing the initial values 
and the updated values set in item 1) of algorithm 1 equals 0. If the result of the determination is affirmative, 
50 item 6) is executed. If the result of the determination is relative, control circuit 11 executes item 2) for calculating 
a new defining point and determining the type, and processing circuit 12 performs the calculation of h2 and 
hi shown in items 3) and 4). respectively, in accordance with the type. If the defining point has been updated, 
control circuit 11 newly selects and updates polynomials of A and G within memory 13 as shown in item 5). 
Control circuit 11 and processing circuit 12 need not be separated, since the above-described processing 
55 may be executed by a software approach by causing the CPU to execute programs corresponding to the re- 
spective control procedures and the above-descrived algorithm. Since the calculation necessary for the above- 
described control and processing comprises multiplication, division, addition and subtraction of simple inte- 
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gers. no particular circuit and processing are required. The above-described operations may be simplified if 
the determination of the type and the flow of the procedure in the control are previously programmed or printed 
in a ROM (read-only memory) and are retrieved (by looking up a table) whenever necessary. 

Since hi and h2 can be simultaneously processed, indepedent processing circuits may be provided for 
5 hi and h2. Furthermore, since hi and h2 are processed in the same manner, M an h2 may be processed by 
a single circuit As will be described later with respect to the effects of the invention, one of the characteristics 
of this algorithm is easiness of performing parallel processing, a plurality of processing circuits may be pro- 
vided. As described above, the circuitry for executing the algorithm 1 can be easily realized. 

10 Second Embodiment 

Next an algorithm equivalent to the algorithm suitable for the Modified decoding algorithm presented by 
Kamiya et al. will be considered. First after the following preparation, the algorithm presented by Kamrya et 
al. wiD be shown. Items which are not described in the following preparation are the same as in the first em- 
15 bodiment 

Preparation 2 (for details, see reference (12)) 

< T : In the present embodiment, the total order is defined as follows: 
Q(i) = a il + b-i2 (a and b are natural numbers which are mutually prime, i = (i1, i2)) 
mSrn, 

if Q(m) < Q(n) is satisfied, or 
only when Q(m) = Q(n) and m2 £ n2. and 
the point next to point m is represented by m ♦ 1 with 
respect to this total order <T. 

2fP(n): = {mel | m2 S n, t s p m, m < r p} 
u: u is a finite partial two-dimensional array having a size q, and is defined as a mapping from Io q (2 (a- 
1)) onto field K. 

F: Bivariate polynomials on field K are expressed by 

f = £fm-z m 

merf, 
where 

z m = x^-y" 12 , rf = {m«=£(a- 1 ) | ftn * 0} 
s = LP(f) = max{m I merf) 
A set of polynomials is expressed by 

f = {fTO. ...fa- 1)}. 

V(u): If the set of un for Q(p) < q is represented by {un I ne7^{2 (a - 1))}, and if df„fl> = 0 (neS.p (a-1+s2)) 
for up, f [up] = 0. At that time, V(up) = {f (a polynomial) | f [up] = 0). 
A : A is obtained in the following manner 

A = UA* ■ » 

q,seIoP(2.(a-1)), 

where A«-» = {mei:(a-1)| mi ^ q1 - s1, m2 =» q2 - s2} when f V(u«) satisfying df q * 0 and LP(f) =« s is present, 

andA*«° <t in other cases. 

h3: h3 of type <i, J> is defined as follows: 

h3 » x rt - »i(D.ffl> - (oVdj)-^ - m * mi - ti(D.g<D 
G: This notation is called a set of auxiliary polynomials of F, and is a set of polynomials: 

G = {g<°>, .-,g<»- D 

Kamiya algorithm 

1) n = (0, 0) , F = {1 , y, y^ . . , y~0 G = «■ 

2) Calculate df n O of all polynomials of F. 

3) If fO satisfying df n fl> * 0 is present a new A and a new defining point t are determined. 

4) The following procedure is executed for all defining points t 
t = (s10>, s20) Provide h3 of type <i, n2 - i>. 
t = (n1 - s1fl + 1, n2 - s20) -► Provide h3 of type <n2 - i, l>. 
t = (n1 + 1, n2-e20>)-> 
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h3 = x" 1 ' » 1 < na - D*i.ftf* I >. 

All polynomials satisfying df„ (D * 0 are removed from F, and all newly obtained h3 are inserted in F. 
5) When A has changed, polynomials of set of auxiliary polynomials G for the new F are selected from 
polynomials of the old G and removed polynomials of F. 
5 6) The process is terminated ifn = n + 1;n = p. The process proceeds to item 2) in other cases. 

An algorithm equivalent to the Kamiya algorithm will now be shown. 

Algorithm 2 



20 



10 1) n = (0. 0). F = {1. y. y-'}. G = 4 „ 

A = {*},B = {u}. 

2) rf the highest-degree coefficients of all polynomials of B does not equal 0, a new A and a new-defining 
point t are determined. 
15 3) The following procedure is executed at all the defining points L 

®t = (S1<0 l s2< , >)-> 

Provide h4 of type (I, n 2 - 1>. 
© t + (n1 - S1C0 + 1, n2 - s20>) 

Provide h4 of type <n 2 - i, i>. 
® t = (n1 +1, n2-s2W)-> 

h4 = X -1{n2-I>.n1-1.fln2-l) > 

All polynomials in which the highest-degree coefficient does not equal 0 are removed from B v and all 
of newly obtained h4 are inserted in B. 

4) Polynomial h3 of the type determined in item 3) is provided. Polynomials corresponding to the poiyno- 
25 mials removed from B are removed from F, and newly obtained h3 is inserted in F. 

5) When A has changed, polynomials of set of auxiliary polynomials G for the new F are selected from 
polynomials of the old G and the polynomials removed from F. Polynomials of set of auxiliary polynomials 
A for the new B are selected from polynomials of the old A and the polynomials removed from B. 

6) The process Is terminated ifn = n + 1;n = p. The process proceeds to step 2) In other cases. 
30 h4: h4 of type <i f j> is def ined as follows: 

h4 = x* 1 " ■ 1 < I >.bG - (Oyd^-X^ " nl * ml - Hffl-aO). 

As in the case of the algorithm 1 and the Sakata algorithm, differences between the algorithm 2 and the 
Kamiya algorithm are in that whether df n O is directly calculated or the highest-order coefficient of polynomial 
b(i) which belongs to B is made to be df n o. Accordingly, if the relationship held between the algorithm 1 and 
35 the Sakata algorithm holds between the algorithm 2 and the Kamiya algorithm, it can be proved that the F ob- 
tained by the algorithm 2 is the same set of polynomials as the F obtained by the Kamiya algorithm in the fol- 
lowing manner. 

If it is assumed thai fflc F of theorem 1 is P> of the Kamiya algorithm, it can be said that the coefficient of 
of the corresponding polynomial W> = flO-u equals df n P) of the Kamiya algorithm. Hence, theorem 1 also 
holds In this embodiment. 
Theorem 2 holds generally. 
Since f q w of theorem 3 equals ft 1 ) of the Kamiya algorithm, theorem 3 also holds. 

Accordingly, as In the case of the Sakata algorithm, since it can be said that w q C0 = b q W f bpW can be up- 
dated by h4. Also in this case, from theorem 1 , d, and dj are the highest-degree coefficients of the polynomials 
45 b n o and b m w, respectively. 

As described above, it has been proved that the F obtained by the algorithm 2 is the same set of polynomials 
as the F obtained by the Kamiya algorithm. 

Accordingly, the algorithm 2 can also be realized by the apparatus shown in FIG. 1. However, the control 
and the processing are different from those of the algorithm 1 . (In the processing, h3 and h4 are calculated. 
50 In the control, the determination of the type is simplified. In memory 13, the polynomials shown in item 1) are 
stored in memory 13 as initial values for the A, B, F and G.) Since the algorithm 2 is also suitable for parallel 
processing, a plurality of processing circu&e may be used. As described above, it can be understood that the 
circuitry for executing the algorithm 2 is easily realized. 

55 Third Embodiment 

Sakata et at. (see reference (1 3)) also show an algorithm for deriving a set of polynomials F for synthesizing 
a multidimensional shift register having a minimum number of storage elements for generating a multidimen- 
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sional array. This algorithm is called a multidimensional BM method. In this embodiment a new algorithm cor- 
responding to the multi-dimensional BM method will be considered. 

Since the detailed algorithm is complicated, characteristics of the algorithm called the BM method will be 
shown in a diagram. The algorithm proposed by Sakata (the multidimensional BM method), which is an exten- 

6 sion of the one-dimensional BM method, including the one-dimensional BM method, has the structure shown 
in FIG. 2. An algorithm corresponding to this algorithm is shown in FIG. 3. As in the above-described two em- 
bodiments, differences between the methods shown In FIGS. 2 and 3 are in that whether df„w is directly cal- 
culated, or the highest-degree coefficient of polynomial bfl> which belongs to B is made to be df n o. Accordingly, 
if theorems 1 through 3 hold also in the multidimensional array, it can be proved that the F obtained by the 

10 multidimensional BM method shown in FIG. 2 is the same set of polynomials as the F obtained by the algorithm 
shown in FIG. 3. 

In the above-described embodiments, points are defined on the two-dimensional space in which n = (n1, 
n2), and each coefficient of the array and the polynomial is a mapping at that point When this cocept is ex- 
tended to the N-dimensional space, points are defined as n = (n 1 . n2 . • • nN). and each coefficient of the array 

is and the polynomial is a mapping at that point Such relationship is identical irrespective of dimension. Hence, 
it is apparent that theorems 1 through 3 also hold in the N-dimensional space. 

Accordingly, it can be said that the F generated from the algorithm shown in FIG. 3 is the same as the F 
generated from the algorithm shown in FIG. 2. Hence, the algorithm shown in FIG. 3 can be easily realized by 
the apparatus shown in FIG. 1. 

20 However, in the case of one dimension, the algorithm shown in FIG. 3 becomes the above-described al- 

gorithm called the Eu method used in decoding RS codes and BCH codes. As will be described in the effects 
of the invention, the algorithm shown in FIG. 3 extends the Eu method to multiple dimensions without loseing 
the characteristics of the one-dimensional Eu method. Accordingly, the algorithm having the configuration of 
FIG. 3 including the algorithms 1 and 2 is termed a multidimensional Eu method which corresponds to the mul- 

25 tidimensional BM method. 

As described above, since the one-dimensional BM method and the one-dimensional Eu method are fre- 
quently used for RS codes and BCH codes, various kinds of modified algorithms have been proposed. A con- 
tinued fraction method (see LR. Welch and R.A.Scholtz: "Continued Fractions and Berlekamp's Algorithm 1 ', 
IEEE Trans. Inf. Theory, IT-25, pp. 19 - 27, Jan. 1979) is a well-known modified algorithm of the one-dimen- 

30 sional Eu method. In this method, the calculation of expression (6) is not performed to the end, and calculation 
not relating to the final result is omitted. However, since this method is the same as the one-dimensional Eu 
method in that B is updated by the calculation of expression (6) and df n <Q is not directly calculated, it is apparent 
that this method is also effective for various kinds of modified algorithms including the continued fraction meth- 
od. 

as As can be understood from the cases of the algorithms 1 and 2, this method does not depend on the method 

of classifying types. Accordingly, this method is also entirely effective for a method in which df n ro Is not directly 
calculated for the multidimensional array, as in the BM method. B is updated by the calculation of expression 
(6), and df n o is calculated as the highest-degree coefficient of polynomials of B. 

Although points and respective parameters (rs. rt and the like) are assumed to be integers, and f„<o and 

40 b n w are assumed to be polynomials, it is apparent that this method is effective even if points and respective 
parameters are extended to arbitrary complex numbers, and polynomials are extended to rational functions 
and the like. 

If the degrees of polynomials u and f are inverted (dual polynomials), the degree of polynomial b = f u is 
also inverted, and d, can be the lowest-degree coefficient of polynomial b. 
45 In theorem 3, it is assumed that d = z v x. However, even if d ■ z** 1 or d = 0, w q <*> « b q <*>. and it is apparent 

that the result of the algorithm having the structure of FIG. 3 does not change. 

Although in the above-described embodiment, the calculation shown in the following expression (7) is used 
for updating B, the same result can be obtained by using the calculation shown in expression (8) except dif- 
ferences in constant terms. 
so - *»-b n <'> - (d/djj.z't.ana (7) 

b n ♦ ,<*> = d r zr«.b n (0 - d,**^©. (8) 

where a„<D = b m u>. 

As described above, differences between the Eu method and the BM method can be described as follows: 
The BM method: df„» is directly calculated, and f^w is obtained using the obtained value. 
55 The Eu method: df n « is not directly calculated, and b„H 00 is obtained as shown in expression (7) using the 

respective highest-degree coefficients d, and dj of polynomials b n ® and b n ®.1^ « is calculated using the same 

d| and d^. 

From such differences, the following differences are produced In providing apparatuses and high speed 
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for the BM method and the Eu method. 

The BM method: In the BM method, it is difficult to calculate both f„« and in parallel. The reason is 

as follows: df n <0 = di is necessary for calculating fn+i 04 . and df n f> is calculated using all the coefficients of f n ro. 

Accordingly, since df ft ro can be obtained only after completing the calculation of f n w. the calculation of f^w 
5 can be started only after the calculation of f„«> has been completed. Accordingly, as shown in FIG. 4(a). in the 

BM method, polynomials f n f> and f^w at points n and n+1 are sequentially calculated. In addition, in the BM 

method, the calculation of f„w. f„^« and df„<o cannot be executed in parallel even if the calculation is performed 

with one dock pulse. Hence, as shown in FIG. 4(b). a calculation time equal to about 2p clock pulses is required. 

This is the characteristic common to the one-dimensional BM method and the multidimensional Bm method 
10 shown in FIG. 2. 

The Eu method: In the Eu method, f n <o and f n *i°° can be calculated in parallel. The reason is as follows: In 
the Eu method, dl is the highest-degree coefficient of b^eB. and polynomials which belong to the B can be 
calculated only using polynomials which belong to B by expression (7). If it is assumed that expression (7) is 
sequentially calculated from the highest degree to lower degrees using b n o> and a^. output b^ 00 is sequentially 

is obtained from the highest-degree coefficient to lower-degree coefficients. Hence, the calculation of polynomial 
bn+2 00 at point n + 2 can be started when the highest-degree coefficient has been obtained (a^® is already 
obtained polynomials). This indicates that polynomials at points n + 1 and n + 2 can be calculated in parallel. 
It is apparent that this also holds for points n and n + 1 . Accordingly, calculation for updating B can be executed 
in parallel at respective points. The calculation of f^**. and f^ <*> can be started when the highest-degree coef- 

20 ficients of bn+1« and b„® have been known. Hence, as shown in FIG. 5(a), the calculation of f^oo can also 
be executed in parallel at respective points. As a result, in the Eu method, if a plurality of processing circuits 
for performing the calculation of expression (7) are provided, a high-speed operation in proportion to the num- 
ber of the circuits can be easily performed. This characteristic is common to the well-known one-dimensional 
Eu method and the multidimensional Eu method shown in FIG. 3, and Is absent in the above-described BM 

25 method. When the calculation of b n <0, b^w, f n (Q and f^w is performed with one clock pulse, as shown in FIG. 
5(b), a calculation time equal to about p clock pulses suffices, and therefore a calculation having a higher speed 
than the BM method can be performed. (In FIGS. 5(a) and 5(b). d, at point n Is represented by d n .) 

In general, a high-speed algorithm indicates an algorithm having a smaller amount of calculation. However, 
for realizing a high-speed operation, an approach to reduce the amount of calculation of the algorithm and an 

30 approach to increase the parallel operations of the algorithm can be considered. This is because the processing 
time can be reduced by simultaneously executing a plurality of processing by parallel processing. A recent pro- 
gress in the VLSI (very large scale integrated circuit) technique can easily provide a chip which can perform 
large-scale parallel processing. Accordingly, in many cases, an algorithm having a high degree of parallel proc- 
essing is suitable for high-speed processing. 

35 As described above, the Eu method of the present invention can perform a higher degree of parallel proc- 

essing than the BM method, and therefore is suitablefor high-speed processing. While the BM method requires 
two kinds of calculation, i.e., the calculation of polynomial f„« and the calculation of scalar quantity df n <0 , only 
the same type of calculation shown in expression (7) is required for f n <» and b„o> In the Eu method. Hence, In 
the Eu method, the processing circuit can be more easily provided in the form of a unit, and therefore it can 

40 be said that the Eu method is more suitable for parallel processing. 

As a result, according to the present Invention, a high-speed processing apparatus which performs parallel 
processing can be easily configured. 

Fourth Embodiment 

45 

Next a description will be provided of an embodiment in which a N-dimensional BM method is used. In 
the N-dimenslonal BM method, addition/subtraction of two points is performed Independently for variables of 
each dimension. For that purpose, in the N-dimensional BM method, more efficient processing can be per- 
formed by using N-dimensional memories each of which deals with variables of the corresponding dimension. 

so Even if calculation relating to (N - 1)-dlmensk>nal variables is simultaneously performed, processing re- 

lating to remaining variables is sequentially performed. Hence, the processes of steps S22 and S24 shown in 
FIG. 2 can be performed in parallel. In order to simultaneously perform calculation relating to (N - ^-dimen- 
sional variables, an Independent memory may be used for each degree of the variables, and respective mem- 
ories may be simultaneously accessed. 

55 Accordingly, in the present embodiment, a method and an apparatus using the method will be proposed, 

in which independent memories for respective polynimlals and respective degrees of variables are realized by 
dividing addresses of a memory or by using a plurality of memories, thereby processing of the multidimensional 
BM method is efficiently performed. 
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FIG. 6 Illustrates an algorithm of the present Invention corresponding to the algorithm shown in FIG. 2. 
FIGS. 7(a) through 7(c) illustrate the operations of the algorithm corresponding to the following fourth through 
sixth embodiments, respectively. In FIG. 6 t the processes of steps S64 and S65 are executed in parallel. 
First the one-dimensional BM method will be considered. In FIG. 8. symbol u represents a memory for 
5 storing a given one-dimensional array u (u t is set for address i). symbol f represents a memory for storing a 
minimal polynomial f (having an initial value of 1) for generating array u. and symbol g represents a memory 
for storing auxiliary polynomial g (having an initial value of 0) for polynomial f. Addresses of these memories 
are controlled from the outside. Since the one-dimensional BM method is considered, only one polynomial is 
stored in each of memories f and g. The (s - 1)-degree coefficient (s is the highest degree of the polynomial) 
10 is set in address I of the memory. Symbol dm represents a register for calculating df n w, symbol df m represents 
a register for holding the value of the df n , symbol IN represents a circuit for outputting the reciprocal of an input 
value, and symbol dd represents a register (having an initial value of Uo) for holding the value of (d^ro/dr,.®). 
Symbol SW represents a switch for selecting an output from memory f only when the defining point has 
changed. Switch SW is controlled from the outside. Symbols x and + represent a multiplier and an adder, re- 
ts spectively. 

At point n, the defining point of polynomial f is represented by s<», and the defining point of polynomial g 
is represented by s«>. Numbers are sequentially added from address 0 in memories f and g, which are controlled 
up to addresses s<o and s®, respectively. At that time, an external circuit compares the values of rs ■•■ sfl> and 
rt + s<», starts to operate a memory having a larger value, and operates the other memory after a time period 
20 corresponding to the number of clock pulses equal to the difference. As described above, since register dd 
holds the value of (dfn (I> /df ra O), respective coefficients of polynomial f at point n calculated and updated by h 
in-step S64 are sequentially input to memory f from address 0 to address sW (s< k > is the updated defining point 
of the polynomial). 

At that time, memory u is controlled such that numbers are sequentially subtracted from address n ♦ 1 
25 until address n + 1 - s<« in accordance with respective coefficients of polynomial f, and respective elements 
of one-dimensional array u stored thereina are sequentially output Thus, df^ 00 at point n + 1 in 2) is calculated 
in parallel with f at point n in step S64. At that time, the preceding df„« is input to register df m . By multiplying 
the reciprocal of the value provided by IN by dW*. value (dfn+i<»/df n <0) can be held in register dd for the sub- 
sequent calculation. However, when the defining point has changed, switch SW is opened, the output of mem- 
30 ory f Is sequentially input to memory g from address 0 to address s«>, and polynomial g is updated. 

Accordingly, in FIG. 8. dock pulses equal to sW is required for one updating operation (one cycle) of poly- 
nomial f. By repeating this operation p times, minimal polynomial f for generating u can be obtained. 

Symbol s^ is the minimum number of clock pulses when df^w and respective coefficients of polynomial 
f are sequentially processed. Hence, in the present invention, useless processing clock pulses are absent in 
35 contrast to the apparatus of reference (11), and therefore the speed of processing is increased. Since high 
speed is not required for the determination and comparison of the defining point, these operations are realized 
by an ordinary CPU. Accordingly, an external circuit can be easily realized by a CPU and a counter circuit for 
performing address control. 

40 Fifth Embodiment 

The external circuit for performing address control of the memories shown in the fourth embodiment only 
omits useless processing clock pulses in the one-dimensional BM method. In a BM method for at least two 
dimensions, however, such an external circuit has an Important role for realizing the above-described multldi- 
45 mensional memories by performing address control and for executing an efficient multidimensional BM meth- 
od. 

First in order to simplify explanation, a two-dimensional BM method wiO be considered. In FIG. 9. symbol 
U represents a memory for storing a given multidimensional array u, symbol F represents a memory for storing 

a plurality of polynomials fw (i = 0 1 - 1) for generating array u, and symbol G represents a memory for 

50 storing a plurality of auxiliary polynomials g© (j ■ 0 1 - 2) for F. Symbol df„f> represents a register for cal- 
culating df„<0, and symbol df m <D represents a register for holding the value of the above-described df„o. A plur- 
ality of such registers are required since a plurality of polynomials are present Other components are the same 
as those shown In the first embodiment 

If it is assumed that i represents the number of polynomials 1», j represents the degree of y*. and k repre- 
ss sents the degree of x*^* (sf> is the defining point of polynomial flO). memories F and G can perform address 
division as shown in FIG. 1 0. and the address can be expressed by (I. j. k). Symbol sj shown In FIG. 1 0 repre- 
sents the highest degree of x for y*. At that time, if it is assumed that the defining point of one polynomial f(i) 
in the set of polynomials F Is sW =» (s1«>, s2W). and a = (a1, a2) (ae I**®) \ 3 a variable, polynomial f» is stored 
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at the position of address (i, s1<Q - a1 , a2). 

If it is assumed that the defining point of one polynomial g© of the set of polynomials G is s© = (s1®. s2<D). 
polynomial g0> is stored in address (j, s1® - a1 . a2) for a e So 8 ® In this case, respective a's used for memories 
F and G need not conform to the total order, and al and a2 can be independently controlled. 

5 For example, first, it may be set such that a1 a a2 = 0, and after operating a1 from 0 to sO, a2 may be 

carried by one, and thereafter a1 may be again operated from 0 to s1. Accordingly. a1 and a2 can be easily 
cent rolled by an ordinary u p-coun ter. As for memory U. address division represented by (J , k) may be performed 
if mapping un of point n = (n1 , n2) is allocated to address (n1, n2), and therefore address control can be easily 
performed as in the case of memories F and G. 

10 Accordingly, the apparatus shown in FIG. 9 may be operated in the following manner First, the external 

circuit compares the magnitude relationship between rs + s<» and rt + s<» with respect to the total order calcu- 
lated from the defining point, operates a memory having a larger value, and outputs polynomials stored in the 
memory. Afterthe time period corresponding to the number of clock pulses equal to the difference, the external 
circuit operates the other memory and outputs polynomials stored in that memory. 

is Since the value of (df n < Vdf m ») is held in register dd, it can be understood that polynomial fw having defining 

point sCO = (s1 w . s2W) calculated and updated by h in step S64 is sequentially input to address (k. s1 (k) -a1 , 
a2)(ae2 0 s <*>) of memory F. At that time, if the point next to point n in the total order is represented by n + 1 = 
n' = (nV. n2') f memory U outputs each element of one-dimensional array u stored in address (nV - a1 , n2' - 
a2)(ae2 0 *< k >) in accordance with variable "a" of updated polynomial fw within the address. 

20 Thus, df n+1 w at point n * 1 in step S62 is calculated in parallel with polynomial fw at point n in step S64. 

At that time, df^ w calculated in step S62 is stored in register df„w, and the preceding df n « is stored in register 
df m . The stored values are used in the succeeding calculation whenever necessary. Therefore, the value of 
(dfnW/dfm©) necessary for the next calculation is held in register dd. and the algorithm shown in FIG. 6 is con- 
tinuously executed. When the defining point has changed, switch SW is opened, the output of memory F at 

25 that time is sequentially input to address (i. siw -a1 , a2)(aeZ^>) of memory G, whereby the set of polynomials 
G is updated. 

Accordingly, in FIG. 9. sw clock pulses are required for updating one polynomial f<D (one cycle), in the multi- 
dimensional BM method, polynomial fw calculated in step S14 and k. i and j of polynomials fro and g® used for 
that calculation are determined in accordance with the type of the defining point. When a plurality of polyno- 

30 mials are stored while performing address division in the above-described manner, k, I and j of the polynomials 
can be sequentially selected by assigning addresses. Accordingly, in order to update all polynomials fw as- 
signed at point n, the number of clock pulses corresponding to the sum of respective s<*> is required. By re- 
peating the above-described operation up to point p, the set of minimal polynomials F for generating two-di- 
mensional array u can be obtained. 

35 It is apparent that the selection of the polynomials can be executed by an external circuit for determining 

the defining point. Accordingly, by selecting a plurality of polynomials by address division using an external 
circuit for controlling addresses, the problem 1) which has been difficult in the conventional apparatus of ref- 
erence (11) can be solved. Furthermore, the problem of the complication of memory access relating to multi- 
dimensional variables described in item 2) is facilitated and can be solved by performing address division for 

40 each degree of y, as shown in FIG. 10. 

Since each polynomial is updated by the minimum number of clock pulses S«, it is apparent that efficient 
processing as in the first embodiment, in which useless processing clock pulses are not provided, is performed. 

It is also apparent that the circuit shown in FIG. 9 is also effective for the N-dimensional BM method if the 
number of division of addresses is increased. 

45 

Sixth Embodiment 

In the fifth embodiment, since a plurality of polynomials are stored while performing address division of 
a single memory, only one polynomial and only one degree can be selected, and parallel processing of poly- 
so nomials and variables shown in problems 3) and 4) is not performed. Accordingly, in the present embodiment, 
an apparatus for realizing parallel processing of polynomials and variables shown in problems 3) and 4) will 
be considered. 

FIG. 11 illustrates an embodiment of the present invention in which parallel processing of polynomials is 
realized. In FIG. 11 , each of memories F and G comprises a plurality of memories for storing a plurality of poly- 
55 nomials in parallel. The thick line shown in FIG. 11 indicates a multiplex signal line corresponding to the number 
of the plurality of memories. Symbols + and x indicate a plurality of adders and multipliers, respectively, pro- 
vided in parallel in accordance with the multiplex signal line. Registers df„«. dfaP. dd and reciprocal generation 
circuit IN are also provided in parallel in accordance with the multiplex signal line. However, sice only one kind 
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of u is stored in memory U and is common to a plurality of polynomials fro, the input line to memory U comprises 
a single line, and the output of memory U is input common to multipliers provided in parallel in common. 

At that time, if addresses of memories F and G are controlled for polynomials having the largest defining 
points at each point n. other polynomials can be controlled by the address control. Accordingly, updating of 

5 all polynomials at point n can be performed by the number of clock pulses equal to s*o, which is the largest 
defining point at that point By repeating the above-described operation p times, the set of minimal polynomials 
F for generating u can be obtained. 

An apparatus for performing parallel processing of variables can also be realized according to the embodi- 
ment shown in FIG. 11. In this case, if divided addresses for respective degrees of y of the polynomial of the 

10 fifth embodiment are stored in parallel in a plurality of memories, calculation relating to variable y within one 
polynomial can be performed in parallel. Hence, parallel processing not only for polynomials but also for vari- 
ables within one polynomial can be realized. 

It is apparent that the circuit shown in FIG. 11 is also effective for the N-dimensional BM method if the 
number of memories is increased. 

is The number of memories of the apparatus of the sixth embodiment is not necessarily equal to the number 

of polynomials or the number of degrees, but a plurality of memories which can be provided may be subjected 
to address division in combination with the case of the second embodiment 

The circuit for performing the processes of steps S62 and S64 can be easily configured by a CPU and the 
like, and is not limited to the circuit configuration shown in the present embodiment 

20 The method of storage, the method of access for memories, and the like which have been described above 

are. of course, efficient not only for the multidimensional BM method, but also for other multidimensional arrays 
and methods for generating multidimensional polynomials. 

In the present invention, since one cycle can be executed with a minimum number of clock pulses s<*> nec- 
essary for sequential processing, useless clock pulses are not provided, and the speed of processing is in- 

25 creased. 

Furthermore, by controlling a plurality of polynomials and variables thereof using address division, the 
present Invention can also be easily applied to the multidimensional BM method. 

In addition, by controlling a plurality of polynomials of the multidimensional BM method in parallel while 
dividing the polynomials in a plurality of memories, the plurality of polynomials can be processed in parallel, 
30 and the processing speed can be increased. 

Moreover, by storing divided addresses for respective degrees of variables within each polynomial of the 
multi-dimensional BM method in a plurality of memories in parallel, processing for respective variables within 
one polynomial can be performed in parallel, and therefore the processing speed can be increased. 

Although the present invention has been described in its preferred form with a certain degree of particu- 
35 larity, many different embodiments of the invention can be made without departing from the spirit and scope 
thereof. It is to be understood that the invention is not limited to the specific embodiments thereof except as 
defined in the appended claims. 



40 Claims 

1. A polynomial-set deriving apparatus for obtaining a set of minimal polynomials for generating a given mul- 
ti-dimensional array, said apparatus comprising: 

first storage means for storing a first set of polynomials to be obtained; 
45 second storage means for storing a first set of aixiliary polynomials for the first set of polynomials; 

third storage means for storing a second set of polynomials different from the first set of polyno- 
mials; 

fourth storage means for storing a second set of auxiliary polynomials for the second set of poly- 
nomials; 

so first discrimination means for discriminating if the coefficient of a predetermined degree of each 

polynomial of the second set of polynomials stored in said third storage means equals zero: 

determination means for newly determinating a defining point if a polynomial, in which the coeffi- 
cient of a predetermined degree does not equal zero, is present as a result of the discrimination by said 
first discrimination means; 

65 first derivation means for deriving polynomials which belong to the second set of polynomials based 

on the value of the defining point determined by said determination means, the second set of polynomials 
stored in said third storage means, and the second set of auxiliary polynomials stored in said fourth stor- 
age means; 
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first updating means for erasing all polynomials, in which the coefficient of a predetermined degree 
does not equal zero, from said third storage means, and for storing the polynomials derived by said first 
derivation means in said third storage means; 

second derivation means for deriving polynomials which belong to the first set of polynomials based 
5 on the value of the defining point determined by said determination means, the first set of polynomials 

stored in said first storage means, the first set of auxiliary polynomials stored in said second storage 
means, and the coefficient of a polynomial of the second set of polynomials relating to the determined 
defining point; 

second updating means for erasing all polynomials corresponding to the polynomials erased by 
10 said first updating means from said first storage means, and for storing the polynomials derived by said 

second derivation means in said first storage means: 

second discrimination means for discriminating the presence of a change of the defining point; 
third updating means for updating the first set of auxiliary polynomials stored in said second storage 
means based on the first set of auxiliary polynomials stored In said second storage means, and the poly- 
15 nomlals erased by said second updating means, when said second dlscrimation means has discrimated 

the presence of a change in the defining point; and 

fourth updating means for updating the first set of auxiliary polynomials stored in said second stor- 
age means based on the second set of auxiliary polynomials stored in said fourth storage means, and 
the polynomials erased by said first updating means, when said second discrimation means has discri- 
20 mated the presence of a change in the defining point 

2. A polynomial-set deriving apparatus according to Claim 1, wherein the updating operation by said first 
updat-ing means and the updating operation by said second updating means are executed in parallel. 

25 3. A polynomial-set deriving apparatus according to Claim 1 , wherein said f irst updating means and said sec- 
ond updating means are realized by common circuit means. 

4. A polynomial-set deriving apparatus according to Claim 1 , wherein said second updating means performs 
the following calculation: 

b q oo = zr»-b n « - (d/dj)-z*b m ©. 
™ where b n ro and b m ® are polynomials which belong to the second set of polynomials, and d, and dj are the 

coefficients of predetermined degrees of the corresponding polynomials b^ and b m ®. respectively. 

A polynomial-set deriving apparatus according to Claim 1 , wherein said second updating means performs 
the following calculation: 

b q W = dj.z™*b n 0> - d r zrt-b m 0), 
where b n o> and b m o> are polynomials which belong to the second set of polynomials, and d, and dj are the 
coefficients of the highest degree or the lowest degree of the corresponding polynomials b n 0> and b m ®, 
respectively. 

6. A method for obtaining a set of minimal polynomials for generating a given multidimensional array, said 
method comprising the steps of. 

setting initial values for a first memory for storing a first set of polynomials to be obtained, a second 
memory for storing a first set of auxiliary polynomials for the first set of polynomials, a third memory for 
storing a second set of polynomials different from the first set of polynomials, and a fourth memory for 
storing a second set of auxiliary polynomials for the second set of polynomials; 

discriminating if the coefficient of a predetermined degree of each polynomial of the second set of 
polynomials stored in the third memory equals zero; 

newly determining a defining point if a polynomial, in which the coefficient of a predetermined de- 
gree does not equal zero, is present as a result of the discrimination; 
80 deriving polynomials which belong to the second set of polynomials based on the value of the de- 

termined defining point, the second set of polynomials stored in the third memory, and the second set of 
auxiliary polynomials stored in the fourth memory; 

erasing all polynomials in which the coefficient of a predetermined degree does not equal zero, 
and updating the third memory by storing the polynomials derived by first derivation means; 
£5 deriving polynomials which belong to the first set of polynomials based on the value of the deter- 

mined-defining point the first set of polynomials stored In the first memory, the first set of auxiliary poly- 
nomials stored in the second memory, and the coefficient of a polynomial of the second set of polynomials 
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relating to the determined defining point 

erasing all polynomials corresponding to the polynomials erased by the updating operation of the 
third memory, and updating the first memory by storing the derived polynomials; 
discriminating the presence of a change of the defining point 
5 updating the first set of auxiliary polynomials stored in the second memory based on the first set 

of auxiliary polynomials stored in the second memory, and the polynomials erased by the updating op- 
eration of the first memory, when It has been discriminated that the defining point has changed by the 
discrimination operation; and 

updating the first set of auxiliary polynomials stored in the second memory based on the second 
10 set of auxiliary polynomials stored in the fourth memory, and the polynomials erased by the updating op- 

eration of the third memory, when it has been discriminated that the defining point has changed. 

7. A polynomial-set deriving method according to Claim 6, wherein the updating operation of the f irst set of 
polynomials and the updating operation of the second set of polynomials are executed in parallel. 
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8. A polynomial-set deriving apparatus for obtaining a set of minimal polynomials for generating a given multi- 
dimensional array, said apparatus comprising: array storage means for storing a given multidimensional 
array u; 

first polynomial storage means for storing a set of polynomials F to be obtained; 
second polynomial storage means for storing a set of auxiliary polynomials G for the set of poly- 
nomials F; 

first calculation means for obtaining polynomials f» which belong to the set of polynomials F based 
on polynomials f» stored in said first polynomial storage means, polynomials g<D stored in said second 
polynomial storage means, and deviations of polynomials df n «; 

second calculation means for obtaining deviations of polynomials df^w based on the coefficients 
of the polynomials fl*> calculated by said first calculation means and the multi-dimensional array u stored 
in said array storage means; and 

control means for controlling accessing operations for said first polynomial storage means and said 
second polynomial storage means and accessed addresses In parallel depending on the degrees of the 
polynomials fw, 

wherein the calculation by said first calculation means and the calculation by said second calcu- 
lation means are executed in parallel. 

9. A polynomial-set deriving apparatus according to Claim 8, wherein said first polynomial storage means 
or said second polynomial storage means includes a plurality of memories, and wherein a plurality of poly- 
nomials are stored in different memories. 



10. A polynomial-set deriving apparatus according to Claim 8, wherein said first polynomial storage means 
or said second polynomial storage means Includes a plurality of memories, and wherein respective vari- 
ables and respective degrees of a plurality of polynomials are stored in different memories. 

40 

11. A polynomial-set deriving apparatus according to Claim 8, wherein said control means controls said first 
polynomial storage means and said second polynomial storage means so that a plurality of polynomials 
are stored in different addresses. f 

45 12. A polynomial-set deriving apparatus according to Claim 8. wherein said control means controls said first 
polynomial storage means or said second polynomial storage means so that respective variables and re- 
spective degrees of a plurality of polynomials are stored in different addresses. 

13. A polynomial-set deriving apparatus according to Claim 8, wherein said first calculation means performs 
50 the following calculation: 

fW = z»-ffl - (df^/df^.z^g®. 

14. A polynomial-set deriving apparatus according to Claim 8, wherein said second calculation means per- 
forms the following calculation: 

df n4 . 1 « = f/U mM 4,. a , 
based on polynomials fW = Z fnW-ZaCaiculated by said first calculation means and the multidimensional 
array u stored in said array storage means. 
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15. A method for obtaining a set of minimal polynomials for generating a given multidimensional array, said 
method comprising the steps of: 

storing a given multidimensional array u in an array memory: 

obtaining polynomials fw which belong to a set of polynomials F to be obtained, based on polyno- 
mials fro stored in a first polynomial memory for storing the set of polynomials F, polynomials g® stored in 
a second polynomial memory for storing a set of polynomials G different from the set of polynomials F, 
and deviations of polynomials dW 

obtaining deviations of polynomials df^w based on the polynomials f*> obtained by the first cal- 
culation operation and the multidimensional array u stored in the array memory; and 

controlling accessing operations for the first polynomial memory, the second polynomial memory 
and accessed addresses in parallel depending on the degrees of the polynomials fw, and executing the 
first and second calculation operations in parallel. 

16. A method of decoding data received from a storage medium or communication channel, the method in- 
cluding deriving a polynomial set using a method or apparatus according to arty preceding claim. 

17. An apparatus for decoding data received from the storage medium or communication channel, the appa- 
ratus including means for deriving a polynomial set in accordance with any of claims 1 to 15. 
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